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"^" ■ Abstract 

f"^ I We explore the possibility of decelerating and Doppler cooling an ensemble of tree-level A-type 

(N ■ 

—^ ' atoms by a coherent train of short, non-overlapping laser pulses. We show that A- atoms can 

jys ■ be Doppler cooled without additional repumping of the population from the intermediate ground 

0^ , state. We derive analytical expression for the scattering force in the quasi-steady-state regime and 

CN . 

analyze its dependence on pulse train parameters. Based on this analysis we propose a method of 

^' 
^1^' choosing pulse train parameters to optimize the cooling process. 
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I. INTRODUCTION 



Doppler cooling [l| relies on radiative force originating from momentum transfer to atoms 
from a laser field and subsequent spontaneous emission in random directions. Cooling by CW 
lasers has been widely studied both theoretically and experimentally within the last several 



decades 



^ . Schemes for cooling the two- level atoms by the trains of ultrashort laser pulses 



-l9| were proposed. The possibility o: 



of counter propagating vr-pulses 
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stimulated cooling the two-level atoms by the pairs 
13| and similar idea of cooling by bichromatic standing 
wave 13| were studied both theoretically and experimentally. The interest in cooling by 
the pulse trains is stimulated by the rapid development of a pulsed laser technology and 



frequency combs (FC) |14|-|16|. In particular the mechanical action of FC on atoms was 
observed experimentally in Ref. |l7| . 

In many cases the atom can not be approximated as a two-level system because the excited 
state may decay to some intermediate sublevels. As an example, group III atoms have no 
single-frequency closed transition on which the cooling of the ground state could be based, 
because their ground states are composed of two fine-structure sublevels nPi/2 and nP-^/2. 
CW laser cooling of this type of A-systems in the presence of bichromatic force-assisted 
velocity-selective coherent population trapping has been studied in [l8|. Other schemes of 
CW sub-Doppler cooling of three-level atoms based on velocity-selective coherent population 
trapping have been proposed earlie r |l9l . |20| . There were also proposals for bichromatic force 



cooling of three-level A-atoms 2l|, |22 |. 



Here we propose the scheme for decelerating and cooling the three-level atoms with the 
ultrafast pulse train. In our scheme both ground states of the A-type system are coupled to 
the excited state by the same laser field. As a result, the cooling does not require additional 
repumping of population from the intermediate state. The exerted scattering force depends 
on atomic velocity via the Doppler shift. Similar to the case of a two-level system, studied in 
[9| , the spectral profile of the scattering force mimics the periodic structure of the frequency 
comb (FC) spectra. Since the positions of FC teeth depend on the pulse-to-pulse carrier 
envelope phase offset, CEPO, the velocity-dependence of the scattering force can be varied 
in time by simply changing the phase offset between subsequent pulses. Thereby, continuous 
compression of velocity distribution in velocity space can be achieved. During the pulse- 
train cooling, continuous velocity distributions gravitate toward a series of sharp peaks 



(typically of the Doppler width) in the velocity space, reflecting the underlying frequency 
comb structure. 

There are several motivations for this work. Wide spectral coverage of FC allows one to 
cool the atoms in a broad range of velocities at the same time. In some cases, FC cooling 
could be used for reducing number of required lasers. Cooling setup based on tunable 
FC can be alternative to Zeeman slowers , whose fields may be detrimental for precision 



measurements 
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23| . The presented analysis is applicable for laser cooling in ion storage rings 
25| where the circulating ions are subjected to chopped laser field. 
We start consideration by deriving analytical expression for the scattering force in the 
quasi-steady-state regime (QSS), based on the expression for the density matrix obtained in 



our previous work [26| . In the quasi-steady-state regime the radiative decay-induced drop in 
the excited state population between two pulses is fully restored by the second pulse. This 
regime is similar to the saturation regime in a classical system of two kicked coupled damped 
oscillators. Based on our analytical expressions, we show that the A-system can be Doppler 
cooled without additional repumping of population from the intermediate ground state. We 
analyze the dependence of the scattering force on the FC parameters. Based on this analysis 
we propose a principle of choosing FC parameters for optimal cooling of ensemble of A-type 
three-level atoms. 

For the pulse-train-driven A-system there are two major qualitative effects: "memory" 
and "pathway-interference" effects. Both effects play an important role in understanding 
of the radiative force exerted by the pulse train on the multilevel system. The system 
retains the memory of the preceding pulse as long as the population of the excited state 
does not completely decay between subsequent pulses. This is satisfied for finite values 
of the product 7T, 7 being the excited state radiative decay rate and T being the pulse 
repetition period. Then the quantum-mechanical amplitudes driven by successive pulses 
interfere and the response of the system reflects the underlying frequency-comb structure of 
the pulse train. If we fix the atomic lifetime and increase the period between the pulses, the 
interference pattern is expected to "wash out" , with a complete loss of memory in the limit 

7T ^ 1. This memory effect is qualitatively identical to the case of the two- level system, 

j 

explored in Ref. [9|. 

The "pathway-interference" effect is unique for multilevel systems. The excited-state 
amplitude arises from simultaneous excitations of the two ground states. The two excita- 



tion pathways interfere. The "pathway-interference" effect is perhaps most dramatic in the 



coherent population trapping (CPT) regime |27H30l| where the "dark" superposition of the 



ground states conspires to interfere destructively, so that there is no population transfer to 
the excited state at all. 

This paper is organized as follows. In Section II we derive analytical expression for the 
scattering force exerted on atoms by the pulse train in a quasi steady-state regime. In Section 
III we study the dependence of scattering force on FC parameters and propose a method 
of their optimization. In Section IV we study the process of cooling of thermal beam of 
three-level A-type atoms by pulse train. We demonstrate that in the optimal cooling regime 
the initial velocity distribution evolves to a comb-like profile with sharp equidistant maxima, 
"velocity comb". The width of each peak is determined by the Doppler temperature limit. 
Finally, the conclusions are drawn in Section V. 

II. ANALYTICAL EXPRESSION FOR A SCATTERING FORCE EXERTED ON 
ATOMS BY DELTA-FUNCTION PULSE TRAIN 

A. Delta-function-like pulse model 

I — I 

As in our previous work [26j] we parametrize the electric field of the pulse train at a fixed 

spatial coordinate as 

E{t)=iEpJ2^ositUct + <^^)git-mT), (1) 

•m 

where i is the polarization vector, Ep is the field amplitude, and $„ is the phase shift. The 
frequency Uc is the carrier frequency of the laser field and g{t) is the shape of the pulses. 
We normalize g{t) so that max \g{t)\ = 1, then Ep has the meaning of the peak amplitude. 
While typically pulses have identical shapes and ^^ = fn4>, one may want to install an 
active optical element at the output of the cavity that could vary the phase and the shape 
of the pulses. 

The A-system, Fig. [T], is composed of the excited state |e) and the ground states \gi), 
\g2) separated by A12; the transition frequencies between the excited and each of the ground 
states are oJegi, ^eg2 correspondingly. The single pulse area corresponding to a transition 
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FIG. 1: (Color online ) Energy levels of A-type system and positions of frequency comb teeth. 
The comb is Doppler shifted in the atomic frame moving with velocity v. 



\9j) -^ |e) is 



Oj = n 



peak 



g{t)dt 



(2) 



where f2^^" = -^(e|D ■ i\g) is the peak Rabi frequency expressed in terms of the dipole 
matrix element. As long as the duration of the pulse Tp is much shorter than the repetition 
time, the atomic system behaves as if it was a subject to a perturbation by a series of delta- 
function-like pulses: ri^j^"" g{t) — t- 6j5{t). In this limit, the only relevant parameter affecting 
the quantum-mechanical time evolution is the effective area of the pulse. The optical Bloch 
equations, in rotating wave approximation, may be written in form: 
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Pee = -IPee ~ ^^ K^ ^ ^T) ^{9jlm [e' 
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Pe,, = -|pe,, + ^E^(^-^^)EV^^'"^*^"'^*"*"^(/^-^^p-/'^.^.)' 
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Sjfl.Pee + 7^J2^(^~ nT)(^,.e^('=^^W-^/*-*")p,^e - ^,e-^('==^W~^^*-*" Ve^y ), (5) 



n=0 



where the detunings Sj = Uc — oOeg^ are the detunings of the carrier frequency from the 
frequencies of transitions \gj) — )■ |e). 



The dynamics of three-level A-type system driven by the coherent train of delta-function 
like pulses has been studied in detail in our previous work |26|. Here we employ analytical 
expression for the density matrix in a quasi-steady-state regime from that work. Although 
general expression for the density matrix was presented there, here we restrict ourself to 
the case most commonly realized. If the energy gap between the two ground states is much 
smaller than the frequency of transition from the ground to excited state, then the ratio of 
decay rates 71/72 is proportional to the ratio of relevant dipole matrix elements in the same 
way as the ratio of pulse areas 9i/92. In this case we can use the following parametrization: 
^1/^2 = 71/72 = tanx- Then the post-pulse excited state population (pgg)^ in QSS regime 
reads 
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i.Ple)r = 8e 2 sin^ — sin^ ttk/D 

D = 8 cos 2x I 4 sin^ h sin^ — cos 271/1; J sin ttk sin {f] + ttk) + 

+ cosTTKCos (f/ + ttk) I 4 COS — (cos27r/i; — 5) + (cos 6 + 3)(3cos27rfi; + 1) — 

— 16 sin^ — sin^ ttk cos(4x) 

—4 cosh I -— I I 4 cos^ — cos 2nK + 2 cos — — cos — 5 J . (6) 

In this formula and below we employ the following notation (see also Fig. [1]) 
(i) The effective single-pulse area 



e = ^jel + 9l (7) 

where 6j are the single-pulse areas for the two transitions \gj) — )■ |e), j = 1, 2. 
(ii) Number of teeth fitting in the energy gap hAi2 between the two ground states 

K = Au/Urep ■ (8) 

Notice that k, generally is not an integer number. When it is integer, the two-photon 
resonance condition is satisfied and the system evolves into the dark state. 

(iii) Doppler shifted phase offset between subsequent pulses 

fi = r]{t)-7]{t + T) = {k,v + 6,)T + <l). (9) 



Here v is the atomic velocity and is the carrier-envelope phase offset between subse- 
quent pulses, i.e., (p = $„ — ^m+i in Eq- dl])- These phase parameter will be used to 
characterize the spectral profile of the scattering force. As shown below the density 
matrix of a system and the scattering force are periodic functions of fj. 

(iv) Residual detunings 6j, j = 1,2, between l^f^) levels and the nearest FC modes in 
the reference frame moving with the atom. In general. Si = [f] + 27ml) /T and 62 = 
{fi + 2'KK, + 27rn2)/T, where integers rij are chosen to renormalize the residual detunings 
to the interval —Urep/2 < 6j < Urepf^- 

Eq. ([6]) gives the value of the excited state population just after the pulse. The time evolution 
between the pulses is described by {mT < t < (m + 1) T)) 

Pleit) = iplelre-''. (lO) 

The dependence on the phase offset f] is the result of interference between the elementary 
responses of a system to subsequent pulses (the persistent "memory" of the system). Par- 
ticularly, when 7T — 7- 00, the excited state completely decays between the pulses and the 
interference factor vanishes (the "memory" is erased), 

( ^ \ 4 sin^ (ttk; 
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At equal pulse areas 9i = 62 and decay rates 71 = 72 (x = ^/4) the equation ([6]) can be 
simplified further 
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D' = I cos {tXk) cos [t] + TTk) I COS^ (vTk) cos'' I -j- I ~ COS — I -|- 

iT\ {_.___, (Q\ , 2/0^ . 2 



cosh ( ^ J (^sin^ (^- J + cos^ (^- ) sin^ {^k) ] ] . (12) 



B. Scattering force 

Now we focus on the evaluation of the cooling force. 



Fz = hkc ^ Im[peg,fiegJ • (13) 



The laser field is present only during the pulse, so effectively we deal with a sum over 
instantaneous forces 

F(t) =PrYl ^A^ - ^T) Imfe-'^'^^^W-^^^-^WVes, (i)] k^, (14) 

where k^ is the unit vector along the direction of the pulse propagation. The change in the 
linear momentum of a particle due to the m-th pulse is Ap^ = lime^o+ JmT- ^if)dt. We 
find 

^^ = {{pZ)r - iPTe)l) kc, (15) 

Pr 

where (p^)^ = Peei^nT + e), {p^Ji = PeeijnT — e) {t <^ e ^ T) are the excited state 

population values just before and just after the pulse. 

This result follows from noticing that Ap^ is an integral of a particular combination 

^ 5(t — mT)6jl\iv[e~'^^^''^^^'^^^^^^'^^^'^ Peg-{t)] over time. This combination enters the r.h.s. of 
i 
Eq. ([3]). Then by integrating Eq. ^\M over time we immediately arrive at Eq. (ITSl) . 

Several insights may be gained from analyzing Eq.( TT5l) . 

(i) Eq. flTSl) simply states that the averaged over a big number of cycles single laser pulse 
fractional momentum kick is equal to a difference of populations before and after the 
pulse. 



(ii) As elucidated earlier for CW laser cooling (see e.g., Ref. [3|) radiative decay plays a 
crucial role in maintaining force directed along the laser beam. In the context of the 
pulse-train cooling, Eq.( TT5|) . radiative decay brings down the excited state population 
in the time-interval between the pulses, thus keeping pre- and post-pulse excited state 
population difference negative; this leads to a net force along the direction of the pulse 
train propagation. 

(iii) In the regime when two PC modes match both transition frequencies between the 
excited and ground states, the system evolves into a "dark" superposition of two 
ground states which is transparent to the pulses. The population of the excited state 
in this case and consequently the scattering force are both zero. 

In the quasi-steady-state regime the value of single pulse fractional momentum kick is 

^ = (PL). X (1 - e--) (16) 

Pr 
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and the average scattering force can be represented as 

F -^ 

^ SC rp • 

In a particular case of equal branching ratios 61 = 62 
scattering force reads (this was obtained using Eq. ( IT2|) ) 



(17) 
1/2, the expression for the 



F = — ^ 

->- SC rp 



Ap hkc sinh7T/2sin^ (y) sin^(7rfi;) 
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D' 



(18) 
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FIG. 2: (Color online ) The dependence of the fractional momentum kick Ap/{pr) on the phase 
offset f] at different values of pulse repetition period T. Solid purple line T = 4 ns, dashed purple 
line T = 50 ns. The parameters of the system are: 7 = 0.05 GHz, = vr/4, k = 0.12. 

In Fig. [2] we plot the fractional momentum kick Ap/p^ as a function of the phase offset 
f]. The radiative force (fractional momentum kick) exerted by the train of coherent pulses 
depends on the atomic velocity via Doppler shift fj = {kv + 6i)T + 0. As velocity is varied 
across the ensemble, the maxima of the force would occur at discrete values of velocities 



(vr {2n — k) — (j))/{kcT), n 



-1,0,1,2,... 



(19) 



In other words the fractional momentum kick (scattering force) spectral profile exhibits the 
periodic structure of the comb (see Fig. [2]). As an example, for T = 5ns and Ac = 600 nm 
carrier wavelength, the force peaks are separated by Vn+i~Vn = 2iT/{kcT) = Xc/T = 120 m/s 
in the velocity space. Depending on the temperature of the ensemble, the comb may have 
several teeth effectively interacting with the ensemble. Notice, however, that if 7T ^ 1 



(Fig. m dashed purple line) the teeth structure of the radiative force washes out and the 
atoms experience radiative force even if their velocities are far away from peaks. In this 
case the power stored in the pulse is delivered to the entire ensemble. This is in a contrast 
with highly-velocity selective CW laser, where the interaction window in the velocity-space 
is typically 1 m/s. 

C. Maximum momentum kick 

The scattering force Eq. f lTTl) is linearly proportional to the post-pulse excited state 
population. Therefore the discussion of the excited state population dependence on FC 



parameters in |26| directly applies to the scattering force too. In Ref. [26| we found that 
the maximum of {p'le)r and correspondingly the maximum of fractional momentum kick for 
the case of equal pulse areas 9i = 62 and decay rates 71 = 72 is reached at optimal residual 
detunings 5i = —62 = —'mod{2'Kn°^^,2Ti)/T and optimal parameter k = k"^* determined by 

K°P* = _arccos(x), (20) 

n 

where x is a root of the following algebraic equation: 



A ...4 oo„ „..u ^^ .;.4 , 1 « _ ^ o ^_ _ _ 



16x^ cos^ 32x cosh -- sin^ - + 16 cos 2x^ 4 cos - + 3 cos(e) + 9=0. (21) 

4 2 4 2 V 2 / 

One can show that for the general case of non-equal decay rates, 71 7^ 72, and pulse areas 
^ = — = tanx 7^ 1 and fixed value of parameter k, the optimal residual detunings are 
determined as 



mod{f]"P\2Tr)/T, \mod{f]°P\2rr)\ < tt 

Si= < , (22) 

{mod{f]"P\ 27r) + 27r)/T, \mod{f]°P\ 27r)| > tt 



mod{f]°P^ + 271 K, 27r)/T, \mod{f]°P^ + 2tik, 27r)| < tt 

{mod{f]°P^ + 2ttk, 2tt) + 2tt)/T, \mod{f]°P^ + 2'kk, 27r) | > tt 

Here 



(23) 
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B 



-.opt _ _ g^jf^loj^ _ Tj-^ _|_ 27r?7,, 



A = cos7r/T;(4cos — -(cos27r/t — 5) + (cos 6 + 3)(3cos27rfi; + 1) — 16sin^ — sin^7r/T;cos(4x)), 



B = 8 cos (2%) (4 sin^ — + sin^ — cos 2tck) sin ttk. 



(24) 



At X = 7r/4 the coefficient B in Eq. flMj) vanishes and f]"^* = — k/2 + 2?™, n = 0, 1... 
After substituting fl2^ into the equation for the density matrix ([6]) one can find the optimal 
value of the parameter k°^*, corresponding to the maximum of the post-pulse excited state 
population and consequently maximum fractional momentum kick. 

The value of the single pulse area B, maximizing the value of the fractional momentum 
kick is equal to vr + 27m, n = 0,1... In Fig. [3] (a,b) we show the dependencies of the QSS 
values of the excited state population {ple)r (^^ ^i = ^2, 7i = 72) and corresponding single 
pulse momentum kick Ap/pr on the effective single pulse area 0. Different curves correspond 
to different values of parameter /i = 7T. The values of {ple)r ^^^ Ap/pr were calculated at 
the optimal value of k, determined by Eq. ( 12T1) for each B and /i = 7T. 




FIG. 3: The dependencies of the quasi-steady-state values of the post-pulse excited state popu- 
lation {ple)r ^^^ single pulse momentum kick Ap/{pr) on effective single pulse area at different 
values of /i = 7T: /i = 10 (dashed pink line), fi = 1/2 (dashd blue line), fi = 1/100 (solid purple 
line) and optimal parameters fj = — vrK"^*, where k°^* is obtained from Eq. (j20p . 



At O = vr the optimal value of the parameter k is equal to 1/2, independently on the 
ratio of individual pulse areas 61/62, 



K, 



opt 
e=7r 
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(25) 



For K, = 1/2 and = vr the excited state population and the fractional momentum kick are: 

iplMQ = n,K=^-) = le^^/Vcosh(7T/2), (26) 

2 7T 

^P/{Pr)e=7T,K=7r,fj=-7T/2 = o t^^h(~^)- (27) 

The spectral resolution of scattering force at k = Kopt vanishes as — >■ tt. 



As it was shown in 26| the maximum post-pulse excited state population in three-level 
A-system (with 61 = 62 = 1/2) 9i = 02 = V^t^) is reached at 7T ^ 1 and is equal to 2/3. 
Consequently the maximum of the fractional momentum kick is also 2/3. 

This result can be generalized to the case of unequal pulse areas 61 7^ 62 and branching 
ratios 71 7^ 72, (in case if 61 7^ im, n = 0, 1). Here at = vr and n = 1/2, the three-level 
A-system, which is initially in the ground state |f7i), eventually reaches the QSS with the 
fractional momentum kick expressed as 

, s N _ ^Pmax ^ 2sin^(2x) . . 

^^''''- Pr (62-&i)cos(2x) + cos(4x)-2- ^ ^ 

If the decay rates and pulse areas are 71 = 7 sin x = ^) 72 = 7 cos x = ^, 6*1 = sin x, 
^2 = cos^ X) the maximum fractional momentum kick f l28|) is equal to 2/3. This limit is 
independent on the value of x (^1 7^ ^rn requires x 7^ vrn/2). 

D. Friction coefficient 

In general, one would be interested in both slowing down the atomic beam and com- 
pressing (i.e., cooling) the velocity distribution. Cooling would occur if there is a negative 
velocity gradient of the radiative force Fg^- One may introduce a friction coefficient /3 by 
expanding the force about some velocity v, corresponding to a certain value of parameter 

f]{v), 

Fsciv + Av) ^ Fsciv) - p{v)Av . (29) 

When the friction coefficient is positive /3 > 0, one observes the compression of velocity 
distribution around v. Negative values of /3 lead to heating of the ensemble. In the limiting 
case when the radiative lifetime is much shorter then the pulse repetition period there is 
no interference between the action of subsequent pulses on a system and consequently no 
velocity dependence of the scattering force. The friction coefficient is thereby /5 = and 
while the ensemble slows down, there is no compression of the velocity distribution. 

12 



The friction coefficient of Eq. (1290 may be directly determined from the analytical expres- 
sion for the force (ITj 



2 . 1 7^ . 2 /^®\ • 2/ B cos {rj + IT k) — A sin {rj + ttk) 



(3 = — 16^fc^sinh — sin — sin (ttk 



2 y ^ ^ {Acos{f] + ttk) + B sin (r/ + ttk) + C) ^ 



C = 8 cosh ^ I (cos (f ) + 1 ) (1 - cos(27rfi:)) + M - cos f ® j j j , (30) 



where coefficients A, B are defined in 

For the case of equal decay rates and pulse areas (x = 7r/4), one has 

/3{x = ^/4) = — - sinh — sin^ — sin^(7rK) cos(7r/i;) sm[f] + ttk) x 

cos^ — cos^ (tck) — cos — I . (31) 

This result depends on the effective pulse area, G, the product /i = 7T and k = Ai2/u)rep- 



In Fig. H] (a,b) we plot the dependence of friction coefficient /3 (at 6*1 = 62) (!3T!) on the 
phase offset f/ at different values of 7T and O at k = /topi, optimally chosen for each pair of 
parameters 7T and G. It acquires the maximum value at fj = fjij, 



where 



Vl3 = - COS I — I - TTK, (32) 



a = cos{'kk) I cos"^ I — j cos^{ttk) — cos I — ) ) , (33) 

b = cosh (^) ("cos^ (^^ sm'iTTK) + sin^ (^]]- (34) 



For the case of non-equal pulse areas x 7^ ^/4) 



r]j3 = — sec ""^ 



^2 + 52 



A (C - D) + V2By/CD - 2 (v42 + 52) - C 



2 



) 



D = V8(A2 + 52) + C2, (35) 

where A, S and C are defined in (^^ [5UI) . 

One can see (Fig. H]) that as the pulse repetition rate grows (smaller 7T), smaller values 
of single pulse area G lead to larger values of the friction coefficient /3. 
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FIG. 4: (Color online) Dependence of the friction coefficient (i/hk1 on phase detuning f] at (a) 
7T = 1/2 and (b) 7r = 1/100. Each panel has three curves with different values of pulse area 0, 
= vr/lO (solid purple line), G = 7r/2 (dashed purple line), = vr (dashed pink line). 

Notice however, that at very small values of 7T ^ 1 the momentum kick per pulse 
becomes smaller and the number of pulses needed to decelerate the atomic beam is increased. 

At very large values of 7T ^ 1, while the friction coefficient /3 vanishes, the momentum 
kick Ap reaches its maximum. If the large value of 7T ^ 1 is due to the low pulse repetition 
rate, then the scattering force Fsc = ^p/T also becomes smaller and the overall cooling time 
is increased. 

For TT-pulse and 61 = 62 (x = 7r/4) the Eq. fl32|) reduces to 



/5. 



hkl 



sinh ^ sin^(7r/i;) cos'^(7rK) sin(?7 + tik) 



(36) 



2 (cos'^(7rK) cos(f/ + ttk) — {cos{2'nK) — 2) cosh {\)Y 
At K = I (chosen in order to maximize the scattering force) the friction coefficient /3^ 
vanishes (similar to the case when 7T ^ 1). One can show that the friction coefficient at 
n = 1/2 and O = tt turns to zero for arbitrary finite ratio of individual pulse areas 61/62 
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and decay rates 71/72 (x 7^ ^/4)- 

E. Finding the optimal cooling regime 

Before discussing criteria for the optimal choice of FC parameters (single pulse area and 
pulse repetition rate ) we analyze the dependence of the scattering force profile on parameters 
7T and O at optimally chosen number of teeth k fitting into the energy gap between the 
two ground states. It is worth noticing, that the optimal value of k°p* is defined with an 
accuracy up the integer number, that is the values n°^^ + n, ra = 0, 1.., where k°^* is defined 
from Eq. fl20l) . are also optimal. Analysis in this section is carried out assuming the equal 
individual pulse areas 61 = O2 and branching ratios 61 = &2 (x = 7r/4). 

In Fig. [5] we study the dependence of the scattering force Fsc on phase offset parameter 
fj at optimally chosen k, Eq. (I2U]) . as the single pulse area G and the parameter /i vary. 

One can see that at small /z = 7T in Fig. [5] the maxima of the scattering force is nearly 
independent on the pulse area 6 as long as 6 > /i. However, as 6 is increased the friction 
coefficient becomes smaller. As an example, at 7T = 1/100 the amplitudes of scattering 
force corresponding to = vr/lO and = 7r/2 are the same, but width of the peaks is 
smaller at = vr/lO. 

At higher values of parameter /x the scattering force saturates at higher values of the 
pulse area 0. But the gradient of the scattering force is decreased. 

At very small pulse areas — )■ the scattering force vanishes (as well as the momentum 
kick Ap) regardless of the parameter 7T. 

To summarize, at lower pulse repetition rates 7T ^ 1 and larger values of pulse area one 
can obtain larger momentum kick and smaller scattering force and compression rate. At 
larger pulse repetition rate (7T <^ 1) and properly chosen one can obtain the maximum of 
compression rate, but smaller momentum kick. In the first case the cooling time is increased 
and the compression of the velocity distribution is slow. In the second case the number of 
pulses needed to decelerate the beam is increased and the scattering force velocity capture 
range is decreased. 

To find the opimal cooling regime one has to compromise between the fast slowing of the 
entire ensemble and its velocity distribution compression rate. In case when the scattering 
force rapidly vanishes in the vicinity of its maxima only the atoms within narrow groups of 
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FIG. 5: The dependence of the scattering force on the Doppler-shifted phase offset fj at the optimal 
value of K = K°P*, chosen according to Eq. (12 ip at fixed value of /i = 7T: panel (a)/i = 1/2, panel 
(b) fi = 1/100 and different values of effective single pulse area 0. Different curves correspond 
to the distinct values of single pulse areas: © = jg (solid purple line), © = f (dashed blue line), 
G = TT (dashed pink line). 

velocity are decelerated. Below we show that this problem can be mitigated. The spectral 
dependence of the scattering force can be varied in time, so that the positions of maxima 
follow the center of velocity distribution of the decelerating ensemble. In this case those 
atoms, which initially were outside of the scattering force velocity capture range and not 
decelerated will be eventually captured by the force profile being moved in the spectral 
domain (e.g. by changing the CEPO (p). 

However, if the initial atomic beam is too fast and (or) the velocity distribution is too 
wide, one can be interested initially in slowing down the ensemble, so that the cooling 
distance will be not too large. In this case one would prefer to have a broad scattering force 
profile (wide velocity capture range). The amplitude of the scattering force has to be large 



16 



enough to mitigate the increase of coohng time and consequently coohng distance. This can 
be reahzed at larger pulse areas 6. For example, at yU ~ 1/2, ~ 7r/2 for the atoms with 
the excited state lifetime r ~ 15 ns at the laser field wavelength A = 589 nm, the velocity 
capture range Avcptr is estimated by Avcptr ~ 20 m/s. The scattering force amplitude is 
quite the same as its maximum value, reached at B = vr. 

At 6 = TT the velocity capture range can be extended up to A/rp, where Tp is the duration 
of pulse. For r^ ~ 1 ps and A = 589 nm the maximum velocity capture range is very broad 
Av^p^^ ~ 5.89 X 10^ m/s. 

If the intial velocity distribution is already narrow and (or) the central velocity value is 
not too high, the priority can be given to the fast velocity distribution compression. 

The optimal set values of pulse area and pulse repetition rate can be chosen based on the 
initial velocity distribution, desired velocity compression rate and the limiting factors such 
as given cooling length and the laser power. 

III. EVOLUTION OF THE VELOCITY DISTRIBUTION 

Now we turn to the dynamics of slowing down and cooling an entire atomic ensemble, 
characterized by some velocity distribution f{v,t) (time-dependence is caused by radiative 
force) . 

A. No-cooling theorem for fixed FC parameters 

Suppose that the positions of FC teeth remain fixed in the frequency domain during the 
deceleration. As the atoms slow down, they come in and out of resonances with different FC 
teeth. The gradient of the scattering force changes its sign (see Fig. Hj) as the Doppler-shifted 
phase (velocity) vary. As a result the sustained cooling can not be realized if the positions 
of FC teeth remain fixed in frequency space during deceleration. This can be demonstrated 
as follows. Suppose the parameters of the frequency comb remain fixed. As a result of 
scattering A^ pulses the atom with initial velocity Vi will be decelerated to the final velocity 
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Vf determined from the implicit equation 

2csc2|csc2 7rfi: / ^^ ^ ,7^/40.2 • 4 ©A 

JSVr = Tf, — kA [Vf — Vi) cosh — COS — sm ttk + sm — + 

fc.Tsinh^f y- y i ^! 2 V 4 A) 

cosTTK I cos^ — cos^ TTK — COS — J {sm{kcTvf + Tin) — sm{kcTvi + ttk)) J . (37) 

where Vr = Pr/M is the recoil velocity. This equation was obtained by integrating Eq. (TTSjl . 
Eq. fl37|) implies that the decrement in velocities would vary across the ensemble. Yet if 
we fix the change of velocity equal to the spacing between the teeth, Vf = Vi — Xc/T, we find 
that the required number of pulses A^o (or time NqT), 

2Ac CSC^ I CSC^ TTK ,lTf 4& . 2f -4© 



No = — r— ^ — cosh ^ ( cos — sin (ttk) + sin —j, (38) 

does not dependent on the initial value Vi. This implies that if we start with a certain 
velocity distribution f{v), the entire distribution is uniformly shifted by —Xc/T every Nq 
pulses: /(f) —^No fiy+Xc/T). Thus, the radiative force exerted by PC with fixed parameters 
does not lead to velocity compression — there is no cooling. 

Notice that the above analysis has neglected variation of intensity across comb teeth. 
Also while there is no compression of the velocity distribution, there is a residual heating 
due to atomic recoil (this arises from treatments beyond our model, see, e.g., Ref. |3|). 

B. Cooling via tuning the FC 

In order to compress the velocity distribution, one has to maintain the positive gradient 
of the absolute value of the scattering force in the vicinity of the center of velocity dis- 
tribution. To attain this condition, the scattering force profile has to follow the center of 
velocity distribution, moving towards the smaller velocities (frequencies) during the process 
of deceleration. In other words, the FC tooth closest to the atomic transition frequency cUegi 
(in the reference frame moving with the center of velocity distribution) has to be somewhat 
red-detuned from u^g^ . Tuning the positions of the FC teeth and consequently the scattering 
force profile can be achieved by tuning the phase of pulses during the cooling process Qj . 

Initially, we start with some velocity distribution /(f,t = 0). To optimize the number 
of cooled atoms, we focus on atoms with velocities grouped around the position of the 
maximum of /(f,t = 0), i.e. the most probable velocity Vmp{t = 0). Radiative force will 
cause both the distribution f{v,t) and the most probable velocity Vmp{t) to evolve in time. 
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To maximize the rate of compression, the friction coefficient needs to be kept at its 
maximum value at Vmp{t)- We may satisfy this requirement by tuning the phase offset 
(j){mT) = $((?7i + l)r) — ^{mT) according to 

(j){t) = -{5 + Kvmpit)) T-f^p, (39) 

where f/^, Eq. ( 135|) . depends only on (time-independent) values of 7T, 9 and Ai2/wrep- As 
Vmpit) becomes smaller due to the radiative force, the offset phase needs to be reduced. 
We may find required pulse-to-pulse increment of the phase offset explicitly 

A0T = 0((m + 1)T) - cl^imT) = -'^Mvp) ■ (40) 

When the phase offset is driven according to (HOl) . there is a dramatic change in time- 
evolution of velocities of individual atoms. As the phase offset is varied over time, the entire 
frequency-comb structure shifts towards lower frequencies. As the teeth sweep through 
the velocity space, atomic v{t) trajectories are "snow-plowed" by teeth, ultimately leading 
to narrow velocity spikes collected on the teeth. This emergence of "velocity comb" was 
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discussed in Ref. [9| for two-level system. Formally, we may separate initial velocities into 
groups 

Vmpit = 0) + (27r(n - 1) - f]^) /k,T < v{t = 0) < Vrap{t = 0) + (27rn - f]p) /kj", n = 0,±l... 

(41) 
The width of each velocity group is equal to the distance between neighboring teeth in 
velocity space, 2-71 /kcT. As a result of "snow-plowing", the n^^ group will be piled up at 
Vnit) = Vmpit) + 27m/ kcT. The final velocity spread of individual velocity groups will be 
limited by the Doppler temperature. To = h'y/2kB- 

To illustrate the train-driven time-evolution for the entire ensemble, we consider a ID 
thermal beam characterized by the initial velocity distribution 

/(^,i = 0) = -^exp(--^). (42) 

The most probable Vmp, average Vave and the r.m.s. Vrms values are expressed in terms oivQ 
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FIG. 6: Time-evolution of velocity distribution for a thermal beam subjected to a coherent train 
of laser pulses. Pulse-to-pulse phase offset of the train is varied linearly in time as prescribed by 
Eq. (|40p . N is the number of pulses, (a) Atomic and pulse train parameters are: jT = 0.25, 
Q = 7r/2. The optimal phase detuning is f] = —1.23. The center of initial velocity distribution is 



Vmp = 500 m/s. 



as 



VSi 



mp — V ^Vq, 



V. 
Vave = \l —Vo, 




(43) 



-"mp 



is the most probable velocity at t = 0. A typical time-evolution of the velocity distribu- 



tion is shown in Fig. |6l Local compression of velocity distribution happens near the points 
Vc{t) + Xcu/T, n = 0, ±1.., where Vc is the time-dependent position of velocity distribution 
center. Clearly, velocity distribution, while initially being continuous, after a certain number 
of pulses develops a comb-like profile. This is the "velocity comb" of sharp peaks separated 
by equal intervals Xc/T in the velocity space. 
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IV. CONCLUSION 

In this paper we studied Doppler cooling of a three-level A-type system driven by a train 
of ultra-short laser pulses. Analytical expression for the scattering force was obtained and 
its dependence on the FC parameters was analyzed. The scattering force Fgc is linearly 
proportional to the quasi-steady-state post-pulse excited state population. Its spectral (ve- 
locity) dependence exhibits periodic pattern mimicking the spectrum of the frequency comb. 
The contrast of the spectral profile of Fsc is a function of the ratio between the excited state 
lifetime and the pulse repetition period, the effective single pulse area and the residual de- 
tunings 5j between the frequencies of individual transitions and nearest FC teeth. In a 
particular case when the pulse repetition period is much longer than the lifetime of the ex- 
cited state, the spectral dependence of the scattering force reflects the broad-band spectral 
profile of a singe pulse. 

The residual detunings 5j can be optimized to maximize the scattering force. At optimally 
chosen detunings the maximum of the scattering force is reached at single pulse area equal 
to TT. However for 7r-pulses the spectral dependence of the scattering force is lost and 
consequently the friction coefficient vanishes. To optimize the cooling process one has to 
compromise between maximizing the scattering force and its velocity capture range and 
maintaining the sufficient gradient of the scattering force (friction coefficient). The spectral 
profile of the scattering force and consequently the friction coefficient can be varied in time 
to follow the moving center of the velocity distribution of decelerating ensemble. This can 
be realized by simply tuning the carrier envelope phase offset. Such manipulation enables 
sustained velocity distribution compression as the atoms slow down. As a result, initially 
smooth velocity distribution of a thermal beam evolves into a series of narrow groups of 
velocities separated by Ac/T, so called "velocity comb". 
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